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Abstract.  Computer  modelling  of  space  configurations  of  ion  clouds  consisting  of  different  m/z  ions 
and  dynamics  of  their  motion  in  the  presence  of  external(dipo!ar  and  quadrupolar)  RF-field  were 
analysed.  Forms  of  ion  clouds  in  ion  traps  were  obtained  by  molecular  modelling.  It  was  shown  that 
stratification  of  different  m/z  ions  appears  as  a result  of  interaction  of  ion  clouds  with  RF  fields  and  with 
each  other.  For  ions  with  high  charge  q=  30  e e|ectron  and  mass  m=30  000  m prolon  (electrosprayed  protein 
ions  in  ion  guides  and  ion  traps)  at  ion  densities  ~108  cm'3  stratification  takes  place  even  at  room 
temperatures.  Coalescence  of  ion  resonance  frequencies  for  light  and  heavy  equally  charged  ions  in  Paul 
ion  trap  was  also  demonstrated. 


INTRODUCTION 

Different  types  of  ion  electromagnetic  traps  (such  as  Paul  ion  trap,  linear  quadrupole 
radio-frequency  trap,  electrodynamic  balance  etc.)  and  quadrupole  and  multipole  ion 
guides  are  currently  in  use  in  different  modern  mass  spectrometers  designed  for 
analyses  of  molecular  mixtures  at  high  throughput  conditions  (proteomics)  when  ion 
densities  and  total  ion  numbers  in  traps  and  guides  are  high  (>106)  and  ion-ion 
interaction  is  essential  to  properties  of  these  devices  even  at  room  temperatures.  In 
case  of  linear  quadrupole  it  was  found  [1,2]  that  at  high  ion  current  in  the  mode  of 
operation  when  ions  are  flying  trough  a quadrupole  and  at  high  charge  density  in  the 
accumulation  mode  (ions  are  accumulated  in  the  quadrupole)  large  discrimination  of 
ions  with  different  mass  to  charge  ratios  (m/z)  takes  place.  Besides,  attempts  of 
selective  ejection  by  dipolar  excitation  of  ions  of  particular  m/z  in  the  presence  of  high 
space  charge  revealed  collective  interaction  of  ion  clouds  of  different  m/z.  Attempts 
were  made  to  explain  some  aspects  of  this  complex  behaviour  of  such  a system  in  two- 
dimensional  zero  temperature  approximation  [3,4 J.  The  segregation  of  ions  of 
different  m/z  occurs  similar  to  segregation  of  immiscible  liquids.  Ions  of  smaller  m/z 
are  concentrated  closer  to  the  centre  of  the  trap.  In  case,  when  supplementary  external 
RF  field  (dipolar  or  quadrupolar)  is  applied  to  ion  ensembles,  a complex  relative 
motion  of  ions  with  different  m/z,  interacting  through  Coulomb  forces,  takes  place.  In 
present  work,  we  approximate  an  ion  cloud  by  finite  number  of  charges.  Each  charge 
interacts  as  a whole  with  a field  of  the  trap,  and  also  interacts  through  Coulomb  forces 
with  other  charges.  The  interaction  of  an  ionic  cloud  with  residual  gas  is  described  as  a 
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friction.  The  ordinary  differential  equations  of  ion  motion  are  solved  by  Euler  and 
fourth  order  Runge-Kutta  methods. 


MODEL 

Two  devices  were  modelled:  Paul  ion  trap  and  linear  quadrupole  trap.  An  infinitely 
long  wire  quadrupole  with  an  infinitely  small  wire  diameter  was  used  as  a model  for 
linear  quadrupole  trap  analysis.  To  trap  ions  inside  such  a quadrupole  we  used  two 
point  charges  of  the  same  sign  as  trapped  ions  positioned  on  quadrupole  axis  at  some 
distance  from  each  other.  This  makes  electric  field  tree  dimensional  without  axial  and 
azimuthal  symmetry.  Analytical  solution  for  the  quadrupolar-dipolar  field  combination 
for  this  electrode  geometry  significantly  simplifies  digital  analyses  of  ion  behaviour. 
We  suppose  that  characteristic  ion  cloud  dimension  is  essentially  less  than  the  distance 
between  electrodes.  We  used  a quasi-potential  approximation  for  description  of 
electromagnetic  forces  in  both  traps.  The  value  of  quasi-potential  (Kapitza-Dehmelt 
effective  potential  energy  [5])  is  proportional  to  the  square  of  radial  (and  in  case  of 
Paul  trap  to  both  radial  and  z-displacement)  ion  displacement. 


3-DIMENSIONAL  STATIC  CONFIGURATIONS 

In  the  frame  of  this  model  we  determined  the  forms  of  ion  clouds  for  the  mixture  of 
different  m/z  ions  at  zero  temperature.  Force  acting  on  i-th  charge  consists  of  the 
electric  field  force  from  all  ion  cloud  charges  q t , each  of  which  is  positioned  at  the 
distance  R from  charge  i 


E a -Ekq  k E i k / IR  ,.kl  S. 
the  force  of  the  field  from  two  end  charges  Q 


E end  =EQ  (R  , -z)  / /R  j -zj3 


(here  z-  ±L  are  co-ordinates  of  end  charges)  and  radial  gradient  of  quasi-potential 


U=  - k (x2+y2)/2. 
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Randomly  chosen  initial  ion  positions  and  velocities  are  relaxing  to  equilibrium 
positions.  Effective  friction  of  ions  in  a buffer  gas  was  chosen  so  that  relaxation  time 
was  some  tens  of  a charge  motion  periods  in  quasi-potential  well.  We  have  found 
different  ion  cloud  configurations  with  characteristic  forms:  disc,  flattened  ellipsoid, 
spherical  ball,  elongated  ellipsoid,  and  thread. This  sequence  of  configurations 
corresponds  to  growth  of  the  ratio  of  transversal  quasi-potential  value  to  longitudinal 
one.  The  length  X of  linear  chain  of  charged  particles  as  function  of  ratio  of  particle 
charge  q to  trapping  charge  Q (placed  at  points  -L,  +L ) can  be  represented 
approximately  by  the  formula 


4(Q/q)‘,2=  (l-x2)/x3'2, 

2X=2  x L is  equilibrium  distance  between  two  charges  q placed  on  one  line  with 
trapping  charges  Q.  The  result  of  calculations  coincides  qualitatively  with  result  of 
experiments  [6,1 1]. 


STRATIFICATION  IN  AN  ION  TRAP 

For  both  linear  quadrupole  and  Paul  traps  pronounced  radial  stratification  was 
demonstrated  with  lower  m/z  ions  dominated  near  quadrupole  axis  and  the  trap  centre, 
similar  to  what  was  observed  experimentally  in  Penning  trap  [7,12]  and  in  linear  Paul 
trap  [11].  For  linear  qudrupole  separation  of  ions  according  to  their  m/z  was  also  found 
in  z-direction  (along  quadrupole  axis).  The  stratification  is  the  result  of  the  balance  of 
the  Kapitza-Dehmelt  effective  potential  force  from  the  RF-field  and  the  Coulomb 
force  from  the  space  charge.  The  dimensionless  motion  equations  for  n charges  have 
a form 


dz/dt  =v2 , 


dvjdt  = -z  + Ez/r3 , 


(1) 


(where  summation  is  taken  over  all  charges  in  trap)  and  analogous  equations  for  x,y 
co-ordinates. 

For  Paul  ion  trap  the  units  of  time  t/,  length  /y  and  velocity  vy  are: 

/,-=  2 Mzo2Q/(Vq), 
z,=2Mzn4tf/V2, 
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Fig.  I . Ion  cloud  structure  as  averaged  over  dimensionless  time  interval  from  t =0  to  1=  10000. 

The  beginning  of  separation  at  kinetic  energy  of  ions  corresponding  to  <v.  2>avcr,1l.l.  =3  (kinetic  energy 
300K).  Clouds  of  higher  and  lower  m/z  are  depicted  separately  with  left  cloud  corresponding  to  high 
m/z  and  right  to  lower  m/z.  a-  transverse  view,  b-  longitudinal  view. 


Fig.2.  Ion  cloud  structure  as  averaged  over  dimensionless  time  interval  from  r=0  to  r=10000. 

Almost  full  separation  at  kinetic  energy  of  ions  corresponding  to  <v;  2>a»c,ai;c  <1  (kinetic  energy  <100K). 
Clouds  of  higher  and  lower  m/z  are  depicted  separately  with  left  cloud  corresponding  to  high  m/z  and 
right  to  lower  m/z.  a-  transverse  view,  b-  longitudinal  view. 
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v/=  2/  // . , where  Mis  ion  mass,  zq  is  the  half  of  a distance  between  end  cup  electrodes 
in  the  ion  trap,  13  is  a frequency  and  V an  amplitude  of  RF  field.  The  average  kinetic 
energy  of  ion  longitudinal  motion  is  Mv2/2  <vz2>  average-  Visual  inspection  on 
computer  monitor  of  ion  cloud  structures  obtained  as  a result  of  equations  (1)  solution 
clearly  shows  separation  of  charges  with  different  m/z  at  <vz2>average  S3  (Fig.l).  Then 
for  Q=2n  0.5  106  Hz,  M=3  10 4 Da  , q = 30  elementary  charges,  zo  = 1 cm,  V=300  V, 
we  have  ti=1.41(fs,  zi=10  ' 2 cm,  v/=0.7  10 3 cm/s  and  kinetic  energy  is 
M V 2 /2  < V_-2>  average  ~ 3 00  ° K. 


For  plasma  density  ~108  cm'3  the  kinetic  energy  is  of  the  order  of  averaged 
electrostatic  energy  per  ion  q2/r  (r  is  average  distance  between  ions),  so  plasma 
coupling  parameter  T=q2/(rT)  ~ 1.  The  same  criterion  can  be  find  in  quite  different 
way.  There  is  a gap  between  inner  droplet  and  outer  envelope  in  the  case  of  two 
component  (two  different  m/z  ions)  system.  These  ensembles  of  ions  could  be 
considered  as  charged  liquids  [12].  The  size  of  the  gap  has  the  order  of  magnitude  of 
the  inner  droplet  radius  R.  If  amplitude  of  thermally  excited  surface  waves  of  the 
droplet 


x = (T/4no),/2  N,/2 

( N ~ n2'3  R2  is  full  number  of  surface  oscillations)  is  less  than  the  gap  between  clouds, 
x<R,  then  the  liquids  would  be  separated,  in  the  opposite  case  they  will  be  mixed. 

We  can  suppose  that 


q~z*e,  e is  electron  charge. 


a ~n  q2, 


Then  liquids  would  be  separated,  if 

nq2  R2>TN. 


So  the  condition  is  nearly  the  same  as 

T=  q2/(4na  T)>1, 


Where  a-(3/(4TC  n))1'3  is  Wigner-Seitz  interparticle  distance.  It  is  possible  that  the 
separation  can  take  place  even  at  room  temperatures,  n=10s  1/cm3,  z = ~ 30  as  was 
described  in  [2] 
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COLLECTIVE  OSCILLATIONS  IN  ION  TRAP 


Charged  ion  cloud  oscillations  - both  longitudinal  (z-oscillations)  in  ion  traps  [8-9]  and 
transversal  in  a linear  quadrupole  traps  [1-4]  were  investigated  experimentally.  The 
amplitude  of  the  oscillations  (parametric  or  forced)  may  be  large  - of  the  order  of  the 
trap  dimensions.  Peculiar  feature  of  oscillations  is  that  some  characteristic  ion 
frequencies,  seen  at  small  ion  densities,  disappear  at  higher  ones.  To  model  resonance 
dipolar  excitation  of  ion  motion  in  a Paul  trap  we  applied  excitation  voltage  to  end 
caps  of  RF-only  Paul  trap.  We  have  found  a strong  influence  of  the  space  charge 
value  on  the  ion  selection  process.  Computer  modelling  revealed  strong  correlation  of 
ion  motion  of  different  m/z,  which  explains  a coalescence  of  ion  resonance  peaks, 
observed  experimentally  [9]  for  isotopes  of  rare  earth  ions.  We  studied  the  coalescence 
of  ion  resonance  frequencies  for  3-D  motion  of  3 light  and  9 heavy  (/h/»»/=1.1) 
equally  charged  ions  in  Paul  ion  trap.  The  light  and  heavy  ion  clouds  are  completely 
separated  after  relaxation  period.  The  configuration  is  flat  because  longitudinal  field  in 
Paul  trap  is  twice  as  large  as  the  radial  field.  After  relaxation  an  external  periodical 
longitudinal  force  was  applied  and  friction  force  was  switched  off.  At  small 
electrostatic  ion  interaction  two  resonances  were  observed  with  frequencies  difference 
about  10%.  As  electrostatic  interaction  increased,  this  difference  goes  to  zero. 


CONCLUSION 

Molecular  modelling  of  multicomponent  system  of  charged  particles  in  a RF-traps 
have  shown  the  possibility  of  different  m/z  ion  cloud  separation  even  at  room 
temperature  in  case  of  multiply  charged  ions  and  high  plasma  densities.  Criterion  of 
separation  obtained  from  the  analyses  of  thermally  excited  surface  waves  of  ion  clouds 
coincides  with  T=1  criterion.  Resonance  peak  coalescence  correlate  with  stratification 
and  takes  place  at  nearly  same  f values. 
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